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ABSTRACT 

We study cosmological simulations of early structure formation, including non-equilibrium 
molecular chemistry, metal pollution from stellar evolution, transition from population III 
(popIII) to population II (popll) star formation, regulated by a given critical metallicity, and 
feedback effects. We perform analyses of the properties of the gas, and use the popIII and 
popll populations as tracers of the metallicity. This allows us to investigate the properties of 
early metal spreading from the different stellar populations and its interplay with pristine 
molecular gas, in terms of the initial mass function and critical metallicity. We find that, 
independently of the details about popIII modeling, after the onset of star formation, regions 
enriched below the critical level are mostly found in isolated environments, while popll star 
formation regions are much more clumped. Typical star forming haloes, at z ~ 15 — 10, 
with masses between ~ 10 7 — 1O 8 M , show average SN driven outflow rates of up to 
~ 10 -4 M /yr in enriched gas, initially leaving the original star formation regions almost 
devoid of metals. The polluted material, which is gravitationally incorporated in over-dense 
environments on timescales of ~ 10 7 yr, is mostly coming from external, nearby star forming 
sites ("gravitational enrichment"). In parallel, the pristine-gas inflow rates are some orders of 
magnitudes larger, between ~ 10 -3 — 10 _1 M /yr. However, thermal feedback from SN 
destroys molecules within the pristine gas hindering its ability to cool and to condense into 
high-density star forming regions. Only the polluted material incorporated via gravitational 
enrichment can continue to cool by atomic fine- structure transitions on short time scales, 
short enough to end the initial popIII regime within less than 10 8 yr. Moreover, the interplay 
between the pristine, cold, inf ailing gas and the ejected, hot, metal-rich gas leads to turbulent 
Reynolds numbers of the order of ~ 10 8 — 10 10 , and contributes to the suppression of pristine 
inflow rates into the densest, star forming regions. 
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1 INTRODUCTION 

One of the most outstanding problems in Astrophysics is the 
study of the birth and evolution of cosmic structures, from early 
times to present days. Our basic understanding relies on the ob- 
servations of a Universe expanding at a present-day rate (Ho) of 
~ 70km/s/Mpc, in which the baryonic content (^ob) is only a 
small fraction (~ 4%). A form of unknown "dark" matter repre- 
sents, instead, the bulk of the matter content (Qom), with roughly 
~ 23% the total energy density of the Universe. The remain- 
ing rsj 70% is accounted for by the contribution of the so-called 
"cosmological constant", A, or "dark energy" (Qq a). Recent de- 
terminations of the cosmological parameters (e.g. iKomatsu et alj 
120101) suggest H = 70 km/s/Mpc (h = 0.70 if expressed in 
units of 100 km/s/Mpc), Q m = 0.272, Q b = 0.0456, and 
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^oa = 0.728. The power spectrum of the primordial perturbations 
follows an almost linear scaling, with an index n = 0.96 and a 
normalization via mass variance within 8 Mpc/h radius a 8 = 0.8. 
As a reference, it is common to define the standard ACDM model 
as the one with the following parameters: Ho = 70 km/s/Mpc, 
Qom = 0.3, Qob = 0.04, Qqa = 0.7, Q oto t = 1.0, a 8 = 0.9, 
n = l. Stars and gal axies are thought to be born in su ch a cos- 
mological setting (e.g. lGunn & Gottll 1 9721 ; IPeeblesll 1 974h . It is be- 
lieved that the first stars have formed from collapse of primor- 
dial, metal-free gas (e.g . iGalli & Pallalll998l: iBromm et aljf l 999 ; 
Abel et all ExXi 12002: lo'Shea & NormaJbOOl lYoshida et all 



20081: Ic en & Riquelme 2008) and they have subsequently polluted 



the surrounding medium by ejecting the metals produced in their 
cores. This process is responsible for the cosmic metal pollution 
and the appearance of the first metals in the Universe. In addition, 
the modalities of the following star formation are altered by the en- 
hanced cooling capabilities of metal-enriched gas. Indeed, in a pri- 
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mordial environment the only relevant coolants are H, He and their 


Cowie & McKee 1977; Brookshaw 1985; Sarazin 1988; Monaghan 


derived molecules JSaslaw & ZiDOvlll967b. like H9. HD or HeH. 


1992; Clearv & Monaghan 


1999; de Avillez & Mac Low 


2002; 


They are able to cool the medium only down to temperatures that 


Klessen & Lin 2003; Jubelgas et al. 2004; Monaghan et al. 


2005; 



result in Jeans-masses ~ 10 — 100 M©, and therefore present-day, 
common ~ M© stars cannot be born. Metal or dust cooling (e.g 



A (7) sta rs cannot be born. Metal or dust cooling (e.g . 
2005: iDwek euii1l2007l:ISchneider & Omukai 2010; 



Omukai et a . ^ 7 _ , . , 

Grassietal.ll2010l : lDwek&Cherchnefjl201lh and turbulence ef- 
fects at the mome nt when the baryons become self-gravitating (e.g. 
IClark et alfeoiol) are thought to be able to lower these limits down 
to solar scales, but within a large mass range. 
Thus, the very first (population III, hereafter popIII) star s are sup- 
posed to be (e.g. lWooslev & WeaverJl995l : ISchaerej2002l) massive, 
and short-lived (up to ~ 10 6 yr), while the following (population 
II, hereafter popll) ones have low mass (solar-like) and can live 
much longer (up to 10 9 — 10 10 yr). Cosmological simulations ad- 
dressing the role of metal enrichment in the history of the Universe 
have been perfo rmed, both by adopting standard recipes for star for- 
mation^ (e g.lKatz & Gunnlll99ll:lcenlll992l:ICen & Ostrikedll992l : 
lKatzlll992l : TKatz et alJll99? TsDringel & Hernquistl l2003h and by 



following the detailed numerical network to account for molecular 

I — I 

evolu t ion and early, catas t rophic run-away co l lapse (e. g. Gnedin 
1998: lYoshida et al l [20031: iBromm et~aDl2003l: IRicotti & Ostrikeil 
200ilWise & Abell2008l : lGreif et al.l2010l : lMaio et alfcOOeSLbOOl 



2010). Rece nt numerical simulatio ns including both popIII and 



popll yields (Torna tore et al. [2007) and, additionally, full molec 
ular evolution dMaio et al.ll2010l) seem to suggest that the transi- 
tion between the two regimes is expected to take place after the 
very first metal-pollution events (dominated mainly by O and C). 
Due to the high metal yields of primordial SN explosions, average 
metallicities reach ~ 10 -3 — 10 _2 Z© on very short time scales 
dMaio etal.ll2oToh . and the subsequent phases of the cosmic star 
formation history are dominated by the population II evolution with 
just minor, residual contr i butions from the population III regime 
(e.g. iTornatore et alll2007l : iMaio et alJ201Ch . 
Also the impacts of primordial supersonic gas bu lk flows 
(iTseliakhovich & Hiratall2oTol : iTseliakhovich et al.ll2010h with re- 
spect to the underlying dark-matter distributions we re recently 
studie d in numerica l simulations for the first time by IMaio et all 
( 2010), followed bv lstacv et all j2oToh . and the conclusion was that 
early structure formation is suppressed, with a cascading effect in- 
ducing a time-delay in gas collapse of some 10 7 yr. 
The detectability of the first stars ha s been widely discussed in the 
literature (for a complete review, see lCiardi & Ferrarall2005h « and it 
is generally believed that next-generation teles copes might be able 
to ob s erve the first galaxies at high redshift (e.g . Weinmann & Lilly 
2005 ; ICrowther et al.l20ld : l Johnson et al.l2010l : lDiikstra & W vithe 



2010). An important point to stress is our ignorance about the 
detailed metal-pollution history at early times, when star form- 
ing haloes are "small" (i.e. with maximum masses of the order 
of 10 5 - 10 6 M at redshift z ~ 20, and - 10 8 - 10 9 M 
at z ~ 10). It is reasonable that such small haloes are not able 
to retain all their gas, given the impact from supernova feedback 
et al.l200ll : ISpringel & Hern auist 2003l: IWhalen et all 



(e.g.lAguirre 1. . . , .. ^ .. . _ , 

2008). Another issue is the behavior of pristine gas in the pres- 
ence of metal-rich gas. Detailed studies of phase mixing are still 
needed to address this question rigorously (see e.g. ISpitzeil l 1962; 



In this paper, we will outline the main features of the collapse of 
H2 -cooled pristine gas and the subsequent pollution from early 
structure formation; in particular, we will address the impact of 
metal-rich gas on pristine, molecular gas and the interplay between 
the popIII and the popll regime. The problems we want to address 
are: is early metal pollution efficient in enriching the surrounding 
gas? How much metal mass is expelled from the star formation 
sites? How much metal mass is trapped in dense, star forming envi- 
ronments? On which timescales? Are metals coming from nearby 
star forming regions or are they re-incorporated by the same halo 
where they were produced? After the first star formation episodes, 
can the pristine, molecular gas inflows restore the popIII regime? 
Which are the preferred popIII and popll star formation sites? How 
do different popIII IMF's influence the picture? 
This work is divided as follows: in Sect. [3 we describe the simu- 
lations used and outline the method adopted, in Sect. [3] we present 
our results, and in Sect. 01 we discuss and conclude. 



2 METHODS 

We use the N-b ody/SPH numerical simulations presented in 
Maio et all (120101). and perfor med by using the entropy-conserving 
code Gadget dSpringell l2005h . The code contains a complete set 
of the most relevant, non-equilibrium, chemical reactions (involv- 
ing H, H + , H~, He, He + , He ++ , D, D + , HD, H 2 , H+, HeH + , 
e~), in order to properly follow the early stages of gas cool- 
ing (e.g. lYoshida et al.ll2003l : IMaio et allEM l200l ExH I2O10L 
and references therein). We also include several metal species 
produced from stellar evolution (C, O, Si, Fe, Mg, S, Ca ), ac- 
cording to two different populations (Torna tore et al.l l2007h : the 
popIII stars, described by a top-heavy IMlQ, and the popll stars, 
described by a common Salpeter IMF. The transition between 
the two regimes is determined by the critical metalli city, Z cr it 
(see also ISchneider et al.1 [20031 : IBromm & Loebll2003h . We con- 
sider several scenarios for the possible popIII IMF, by adopt- 
ing different IMF mass ranges, slopes, yields, and critical metal- 
licities. In the cases where a top-heavy IMF was assumed (see 
also e .g. IRicotti & Ostrikedl2004l : lKarlsson et al.ll2008l : IMaio et al.1 
2010), the pollution was faster and stronger, reaching ~ 10 — 
10 -2 Zq immediatelly after the onset of star formation. The cases 
of moderate- or low-mass IMF show a similar, but more gradual, 
trend. Since the goal of this paper is an investigation of the dy- 
namical properties of the pollution process, without loss of gen- 
erality, we will focus mainly on the top-heavy IMF runs, and we 
will only refer to the other runs as comparison. We consider dif- 
ferent simulations in the framework of the cosmological standard 
model, assuming Z cr u = 10~ 3 , 10~ 4 , 10~ 5 , 10~ 6 Z©, using high- 
resolution boxes of 1 Mpc a side, and 2 x 320 3 particles for gas and 
dark-matter. 

The main processes which cause pristine gas to cool are particle 
collisions and molecule creation, with typical molecule fractions 



Several studies use a different way to deal with star formation, i.e. 
via "sink" partic les, but we warn the reader that very recent calculations 
( Bate et al. 2010) have shown two main limitations of this approach: over- 
prediction of the accretion rates, and mispredictions of the particle orbits. 



2 The exact slope for the popIII IMF has little effect on the metallicity 
evolution. The main contribution to the metallicity is via pair-instability 
supernovae (PISN) whose fraction only show slight changes between ~ 
0.35 — 0.40, for variations of the slope within the limits of ~[-3.0, -0.5], 
given an IMF range of ~ [100 M , 500 M ] jMaio et al.ll2010h . 
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becoming larger than ~ 10 -2 and completely dominating early gas 
cooling. At the bottom of the cooling branch, we assume that star 
formation and fee dback processes occur fe.g. lspringel & Hernquistl 
l2003l :l Maio et al. 2009), and metals are produced from stellar evo- 
lution. The feedback processes we deal with are thermal feedback, 
mechanical feedback (both shocks and winds) and the chemical 
feedback: thermal feedback from supernova? injects entropy into 
the medium heating up closeby particles and destroying molecules; 
kinetic wind feedback is stochastica lly taken into account, accord- 
ing to Springel & HernquisJ (120031) . and it is responsible for re- 
moving away from the collapsing sites and spreading around gas 
and metals at velocities of the order of some ~ 10 2 km/s; chem- 
ical feedback rules the actual stellar IMF and leads the transition 
from the primordial popIII regime to the following popll one. In the 
popIII star formation regime, the only stars contributing to metal 
pollution are the ones in the mass range [140 M©, 260 M©], that 
explode as PISN after a life-time of at most ~ 10 6 yr. Their con- 
tribution to the to tal SFR rapidly drops down ~ 10 -3 — 10 -4 (see 
Maio et al"]|201(t) . and this seems consistent with what can be in- 
ferred from more recent observational analyses (e.g. lPapovich et al.1 
l20ld and private communica tion). Given the many uncertainties 
on the primordial IMF (e.g. see Yoshida 2006; Yos hida et al J2Q07I : 
ICampbell & Lattanzioll200 8: Suda & Fuiimotdl201(t) . we also con- 
sider the other limiting case of a Salpeter-like popIII IMF, with 
mass range similar to popll, but with metal yields consistent with 
those from metal-free stars (see Sect. 4.1. in lMaio et alj|201(j for 
further details). 

In the case of popll, the surrounding environment is polluted via 
SNII and SNIa, which occur - 10 7 - 10 8 yr, and - 10 9 - 10 10 yr 
following the formation of their progenitor stars, respectively. 
Moreover, stars can loose a significant fraction of mass during their 
AGB phase. Metals are injected into the ISM and pollute it by 
means of stellar winds. First star formation events take place at 
z rsj 16.3, when the Universe is roughly 2.4 x 10 8 years old. 
In the following, results from our analyses will be presented and 
discussed. We will refer to enriched gas particles with metallicity, 
Z, higher than or equal to Z cr it as popll particles, while to the 
ones with metallicity below Z cr it as popIII particles. Particles with 
Z = are referred to as pristine particles. 



3 RESULTS 

In this section we present our results, paying attention first to some 
global properties of metal enrichment (Sect. l3.lt , then to the statis- 
tical (Sect. l3.2t and dynamical (Sect. l3.3t ones. We also discuss the 
features and the interplay between the two different stellar popula- 
tions (popIII and popll). 



3.1 Volume filling factor 

To have a general idea of the metal spreading process, we start pre- 
senting the basic features of the volume filling factor. 
In Fig. Q] (left panel), the volume filling factor of metal enriched 
particles as a function of the redshift is shown. The volume filling 
factor, fv, is defined as follow^: 



fv 



J2i m z,i/Pz,i Y^i m Z,i/ 'pz,i 



V 



(1) 



with i and j integers running over the number of the enriched par- 
ticles and of all the SPH particles, respectively, mz particle metal 
mass, pz metal-mass density, m particle total mass, p total-mass 
density, and V simulation box volume. The solid line refers to the 
filling factor computed by considering only the polluted particles 
(Z > 0), the dashed line refers to the particles with metallic- 
ity < Z < Z cr it (popIII regime), and the dotted line to the 
particles with metallicity Z ^ Z cr it (popll regime). In all the 
three cases, the general trend is similar: the filling factor increases 
as redshift decreases. This is simply due to ongoing metal pollu- 
tion events at z < 16. In more detail, at early times, some met- 
als are produced, after the first PISN explosions, and the produc- 
tion sites are the very dense, clumped regions (with overdensities 
6 > 10 4 ), so the resulting volume filling factor is rather small 
(rsj 10 -16 ). At later times, particles are ejected from these dense 
regions and their densities rapidly decrease, which results in an in- 
crease of fv- In this regime, low-density particles, which experi- 
enced mechanical (winds) and chemical (metal pollution) feedback 
from stellar evolution, are dominating. Interestingly, the discrimi- 
nation by Z cr it shows a higher filling factor for the particles with 
< Z < Zcrit, at redshift ^ 14 — 16, and a behaviour dominated 
by the Z ^ Z cr %t - particles, at z < 14. This has to be inter- 
preted as a consequence of the metal-enrichment process, as most 
of the enriched regions enter the popll regime very rapidly, while 
the residual popIII areas are confined to isolated regions or on the 
border of popll star forming sites and contribute less to fv- The 
exact value of the critical metallicity, Z cr it, does not alter signifi- 
cantly the properties of metal pollution, because of the high metal 
yields of early stars: larger Z cr it slightly postpone the epoch when 
popll becomes dominant, while smaller Z cr it anticipate it. In the 
extreme case of Z cr %t — 10 ~ 6 Z©, then, fv is immediatelly led 
by the popll regime, basically since the onset of star formation. A 
Salpeter-like popIII IMF will only induce a short delay i n the over- 
all be haviour, due to the longer stellar lifetimes (see also Mai o et all 
l2010i) . To illustrate the insensitivity of the volume filling factor to 
the adopted value of Z cr it, in the same Fig.Q](right panel) we com- 
pare the filling factors obtained for simulations having different 
critical metallicities: dotted, solid, dashed, dot-dashed lines refer 



tO Zcrit 
Zcrit — 



io _d z , z t 



io- 4 z©, z t 



10" 5 Z©, and 



10 Zq, respectively. As expected, there are no strong 



3 Here, there is no difference by defining fy via physical or comoving 
quantities. 



variations among these cases, and despite some differences at the 
initial times, the overall trends easily converge. The one-order of 
magnitude difference in the very early epoch is a consequence of 
the chemical feedback: higher values of Z cr it allow the popIII 
regime to last longer and to pollute more with respect to the lower- 
Z crit cases. 



3.2 Statistics 

To go deeper into the analysis, we study in detail the relevant sta- 
tistical properties of the gas during the early stages of metal enrich- 
ment, and then we show the implications for the popIII and popll 
regimes. 



3.2.1 Global features 

In Fig. |2 we show the redshift evolution of metal-enriched particles 
in the phase diagram (p — T), at z = 15.00 (left), z = 13.01 (cen- 
ter), and z = 11.00 (right). The color scale denotes the fraction of 
particles at any given density and temperature. The corresponding 
age of the Universe is 2.66 x 10 8 yr, 3.24 x 10 8 yr, and 4.09 x 10 8 yr. 
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Figure 1. Left panel: metal filling factors for different regimes. Solid, dashed, and dotted lines are the filling factors computed for particles having metallicity 
Z > 0, < Z < Z cr it, and Z ^ Z cr n (Z cr i t = 1O _4 Z ), respectively. Right panel: metal (Z > 0) filling factors for simulations with different critical 
metallicity. Solid, dotted, dashed, dot-dashed lines refer to Z cr n = 1O _3 Z , Z cr n = 1O _4 Z , Z cr n = 1O _5 Z , and Z cr n = 1O _6 Z , respectively. 
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Figure 2. Redshift evolution of metal-enriched particles in the p — T diagram. Data refer to redshift z = 15.00 (left), z = 13.01 (center), and z = 11.00 
(right). The color scale corresponds to the fraction of particles in a given pixel (normalized to 1). The critical density at z = 11, in the standard ACDM 
model, corresponds to ~ 10 -26 /i 2 g/cm 3 . In the right panel, mechanical and thermodynamical evolution of the gas is shown at different stages, via different 
arrows (see labels): low-temperature cooling from molecules (mainly H2) or fine- structure (FS) transitions leads to star formation (SF) and subsequent gas 
outflows into the low-density regions; there hot gas cools down to ~ 10 4 K via resonant atomic transitions and below 10 4 K via molecular and fine- structure 
transitions. 
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Figure 3. Redshift evolution of gas density distribution. Data refers to redshift z = 15.00 (left), z = 13.01 (center), and z = 11.00 (right). The histograms 
filled by 45° blue lines refer to pristine gas, while the ones by 135° red lines to enriched gas. Normalization is fixed to 1 for the maximum, in both cases. 
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Figure 4. Normalized joint probability distribution (left column), P(S, Z), between metallicity, Z, and over-density, S, at redshift z = 15.00 (top), z = 13.01 
(center), and z = 11.00 (bottom). The corresponding conditional probability distributions, P(Z\S) and P(5\Z), are also shown in the central and right column, 
respectively. In each panel, the horizontal, solid line is drawn in correspondence of the critical metallicity, Z cr u = 10~ 4 Zq and the vertical, solid line marks 
the star formation over-density threshold. 



In the plots, we immediately recognize the gas which falls in the 
haloes and gets shock heated (low densities and low temperatures), 
the isothermal phase (where the heating is couterbalanced by cool- 
ing, resulting in a roughly constant temperature of ~ 10 4 K), and 
the catastrophic cooling branch of over-dense gas (high densities 
and low temperatures) that enters the star forming regime as soon 
as it reaches a threshold of ~ 140 h 2 cm -3 . Above this threshold, 
stochastic star formation and feedback effects take place (high den- 
sity and high temperatures) and lead the hot outflows from dense 
environments into lower-density regions (low densities and high 



temperatures). What emerges from the diagrams is also the be- 
haviour of metal enrichment after the initial pollution events. At 
the beginning (left panel), only the gas at high number densities 
(~ 10 2 - 10 3 h 2 cm ) is affected, as this is where star forma- 
tion is occurring. As star formation and metal spreading proceed, 
the surrounding, lower-density regions are the next to be affected 
by metal pollution. The high temperature particles (T > 10 4 K in 
the diagrams) are the ones impacted by SN feedback: in particular, 
hot, dense gas has been heated, within star forming regions, by the 
high SN temperatures (thermal feedback), while hot, low-density 
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gas has been ejected from star forming regions by SN winds (wind 
feedback). Interestingly, the enrichment pattern of the phase dia- 
gram follows the molecular catastrophic-cooling branch, meaning 
that an "inside-out" pollution mode is taking place. Indeed, since 
the critical density, in the ACDM model, at z ~ 11, is roughly 
10 -26 /i 2 g/cm 3 , Fig. [2 suggests that the very early enrichment 
process affects the collapsed object^ by moving from the center 
towards the borders and the lower-density environments. 
Correspondingly, in Fig. [3] we show the gas distribution (his- 
tograms) at redshift z — 15.00 (left panel), z — 13.01 (central 
panel), and z — 11.00 (right panel). The histograms filled by red 
135° lines refer to enriched gas, and the ones with blue 45° lines to 
pristine gas. Most of the particles at low densities are pristine, while 
the peak of the distribution for metal-enriched gas is at high den- 
sity, where most of the metals are being produced: these "hybrid" 
particles are the ones where the sub-grid modeling of the multi- 
phase medium and star formation is being applied. As time passes, 
metal spreading produces a split, with a second peak at lower and 
intermediate densities. This second peak represents the isothermal 
stage before the run-away collapse. Its location is very meaningful, 
since it gives us some idea of the environment affected by metal 
pollution. The reason for the aggregation of enriched particles in 
this stage of the phase diagram is due to the long cooling time- 
scales at low density, in comparison to the cooling branch. Particles 
that underwent star formation process are spread around from wind 
feedback. They are hot (T ~ 10 5 — 10 6 K) and are ejected into 
lower-density regions, where their cooling capabilities allow the 
temperature to rapidly drop down to T ~ 10 4 K. Subsequent cool- 
ing is slightly slower (as observable from the larger spread below 
~ 10 4 K, in Fig. |2}, due to the less efficient metal fine-structure 
transitions, but it still dominates the total cooling function, with 
respect to the molecular contribution. The lack of particles in cor- 
respondence of high/intermediate densities and high temperatures 
is due to the feedback effects combined with strong cooling, which 
allow the gas to easily cool down to T < 10 4 K. Indeed, due to 
feedback from SN, entropy is injected into the star forming medium 
and the gas is pushed away towards low-density areas. The typical 
velocity for stellar winds is 500 km/s, thus in ~ 10 6 — 10 7 yr the 
gas can be pushed several kpc awa)|j , beyond the isothermal peak 
(see Fig. |3, where densities are orders of magnitude smaller. 
By comparing the distribution at different redshifts, it turns out that 
the mass fraction of enriched material located in over-dense (i.e. se- 
lected with over-density 8 = p/ p cr (z) > 10 2 ) regions is still more 
than rsj 99% at z ~ 13, and decreases down to ~ 87% at z ~ 11, 
i.e. At rsj 0.08 Gyr later. At this redshift, the fraction of new met- 
als, fz,new, represents more than 90%, and the corresponding en- 
richment rate, sz = dfz/dt ~ fz,new/At, is ~ 10 _8 yr _1 . 
When looking at the amount of metals moving along the cooling 
branct0 we find that ~ 1%, ~ 6%, and ~ 4% of the enriched mass 
at redshift z = 15.00, z = 13.01, and z = 11.00, respectively, is 
undergoing run-away collapse. 

Comparison of the simulations with different Z cr it shows that this 
picture does not change dramatically, and the enriched mass un- 
dergoing run-away cooling increases, in all the cases, up to around 

4 According to the "top-hat" model, collapsed objects are the ones with 
densities of the order of ~ 200 times the critical one. 

5 The physical distance travelled at 500 km /s in 10 6 — 10 7 yr is ~ 0.5 — 

5 kpc, corresponding to ~ 8 — 80 kpc comoving, at redshift 15. 

6 The cooling branch roughly corresponds to densities higher than ~ 10 2 
times the critical one and below the star formation threshold of about 
140 h 2 cm -3 ~ 10~ 22 h 2 g/cm 3 . 



5% - 7% at z = 13.01, and reaches about 3% - 6% at z = 11.00. 
If we look at the run with a Salpeter-like popIII IMF (see details in 
Mai o et alJl2O10L in particular Fig. 7), stellar lifetimes are longer, 
and consequently also the time-delay for metal spreading is longer 
( ~ 10 7 yr). So, metal-rich particles undergoing run-away cool- 
ing are significantly present only after z < 14 and account for 
~ 3% — 8% of the enriched mass. 

In Fig. |4] we investigate how gas mass metallicity correlates with 
over-density. We consider metal-enriched particles and compute 
the joint probability distribution between Z and 8, P(8, Z), and 
the conditional probability distributions, P(8\Z), and P(Z\S). The 
joint probability distribution, P(8, Z), is computed by summing 
up all the metal contributions at any 8 and Z and normalized over 
the whole {5, Z} space. While the conditional probability P(Z\S) 
is computed by summing up all the contributions for a given 8, 
and normalized over the corresponding {Z} subspace. Similarly, 
the conditional probability P(S\Z) is computed by summing up 
all the contributions for a given Z, and normalized over the cor- 
responding {8} subspace. The horizontal line marks the critical 
value, Z cr it, while the vertical line marks the star formation den- 
sity threshold (~ 140/i 2 cm -3 ), above which the star formation 
sub-grid modeling is applied. The plots refer to redshift z = 11 
(upper row), z = 13 (central row), and z = 15 (bottom row). 
The joint probability distribution, P(S, Z), has a maximum around 
Z ~ 10 -6 — 10 ~ 4 , almost despite of the redshift considered, and its 
time evolution shows clearly how metals are ejected by wind feed- 
back from the star forming regions, at 8 ~ 10 4 — 10 5 , to the lower- 
density ones, at 8 < 10 2 . This process is better highlighted in the 
conditional probability distribution P(Z\8): at earlier times, during 
the very first episodes of star formation, there are still few metals, 
mostly at low Z, and they have not reached yet the low-density 
environments. Later on, P(Z\S) is rapidly boosted to higher metal- 
licities, with Z overcoming Z cr it at any 8, as a consequence of 
the high metal yields of primordial stars. By checking the condi- 
tional probability P(8\Z), it emerges that, for any fixed Z besides 
the star-forming particles at 8 ~ 10 4 — 10 5 , the regions where 
enrichment due to metal- spreading is dominant are soon the ones 
with 8 < 10 2 (see also previous sections). Furthermore, the over- 
density range between ~ 10 2 and ~ 10 4 , at z ~ 11, presents a gap 
in the distribution, because it corresponds to the fast, catastrophic 
run-away cooling, through which enriched gas rapidly re-collapses 
to higher densities and lower temperatures. 



3. 2. 2 PopIII and popll 

After our investigation of the global features of the primordial en- 
richment, we want to study the differences between population III 
and population II regimes. 

In Fig. [5] we select particles according to their metallicity and plot 
the density distribution of metal-rich particles for both popll and 
popIII ones (in a similar way to Fig.|3}. By comparing the distribu- 
tions of particles that have Z ^ Z cr it (first row) and Z < Z cr it 
(second row), it is easily shown that, in both cases, enriched star 
forming particles are continuously produced by the ongoing star 
formation, and, while the Universe evolves, wind f eedback pushes 
them towards external, lower density regions (see [Tornatore et al.l 
l2007h . In particular, the relative fraction of popIII particles ouside 
the collapsed objects is more than a factor of 2 higher than the popll 
one. Let us discuss these two cases more in detail (see also Table 

ID. 
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Figure 5. Redshift evolution of gas density distributions of pristine and metal-rich particles, selected by their Z. In both rows, we show the pristine-gas 
distribution with histograms filled by 45° blue lines, while enriched gas distributions, for particles having Z ^ Z cr u in the upper row and Z < Z cr u in the 
lower row, are plotted with hystograms filled by 135° red lines. Data refer to z = 15.00 (left), z = 13.01 (center), and z = 11.00 (right). Normalization is 
fixed to 1 for the maximum. 



Table 1. Mass fractions of enriched gas. 



Table 2. Mass fractions of the enriched cooling gas. 



Redshift 


popll 
relative fraction 


popIII 
relative fraction 


Redshift 


popll cooling 
relative fraction 


popIII cooling 
relative fraction 


15 


43% 


57% 


15 


3% 


97% 


13 


94% 


6% 


13 


80% 


20% 


11 


99% 


1% 


11 


94% 


6% 



• Popll. The mass fraction of metals with Z ^ Z cr %t is initially 
zero, but reaches already ~ 43% at z ~ 15 and rapidly rises up to 
~ 94%-99%atz ~ 13-11, and above at later times. The polluted 
material is pushed away from the over-dense regions and expands 
with time. In particular, about ~ 23% of formed popll material has 
travelled into surrounding 6 < 10 2 regions, by redshift z ~ 11 - a 
time lag of ~ 2 x 10 8 yr. 

• PopIII. For the very-low metallicity particles, the picture is 
quantitatively quite different and particles with Z < Z cr it are 
found outside the collapsing regions, in a more significant way. 
The onset of star formation is initially dominated by popIII, but 
shortly after, at z ~ 15, the residual enriched mass fraction with 
Z < Z C rit drops down ~ 57%, and becomes just ~ 6% at z ~ 13, 
and~ l%atz~ 11. 

Let us now focus on the properties of the enriched cooling 
mass only (see Table 13, i.e. the mass which has experienced wind 
feedback and is cooling back into star forming sites. In the initial 
stages, as mentioned in Sect. 13.21 roughly ~ 1% of the enriched 
mass is in the catastrophic cooling regime, and this fraction slightly 
increases up to 4% — 6% at z ~ 11. At z ~ 15, the enriched cool- 
ing mass is made of popIII type for ~ 97% and of popll for ~ 3%. 
After that, the situation completely changes and, at z ~ 13, the 
popIII part is only ~ 20%, and the popll rapidly increases up to 



~ 80%. Later on, at z ~ 11, the residual popIII mass fraction 
drops down to ~ 6% and the popll one increases up to ~ 94%. At 
lower redshifts, the popIII part becomes negligeble with respect to 
the popll one. Also in this case, different values for Z cr it do not 
induce drastic changes: the popIII material is significantly present 
only at the very first stages and completely replaced by popll at 
z < 14 — 13, although occurs more quickly for lower Z cr it. 
When a Salpeter-like popIII IMF is adopted, the residual collaps- 
ing popIII material at z < 14 is again only a few percent (more 
precisely, ~ 3% at z ~ 11). 



3.3 Dynamics 

To get a better understanding of metal pollution at early times, 
we discuss the dynamical features of the material involved in the 
spreading process. As an example, we focus on the evolution of the 
most massive dark-matter halo (mass between ~ 10 7 — 10 8 M / h 
at z ~ 16 — 10) and compute radially averaged quantities. The co- 
moving virial radius is computed as the distance from the center of 
mass at which the average enclosed over-density is ~ 200, and, at 
these stages, it varies between ~ 4 — 9kpc//i. We note that such 
over-densities include the whole cooling branch and a relevant part 
of the isothermal gas phase (see e.g. Fig. O. 
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Figure 6. Mass profiles as a function of comoving radius for enriched (left panel) and pristine (right panel) gas at redshift z = 11.00 (solid lines), z = 13.01 
(dotted lines), and z = 15.00 (dashed lines). Each vertical line corresponds to the virial radius at each redshift. The comoving gas softening length is 
~ 0.1 kpc/h and the minimum gas smoothing length is~0.14kpc//i. 
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Figure 7. Mass outflow rates for enriched gas (solid lines), and inflow rates for both pristine gas (dotted lines) and enriched gas (dashed lines). Redshifts are 
z = 16.01 (upper left), z = 15.00 (upper center), z = 14.02 (upper right), z = 13.010 (lower left), z = 12.00 (lower center), and z = 11.00 (lower right). 
The vertical dot-dashed lines are the corresponding virial radii. 
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Figure 8. Redshift evolution of the velocity field lines for pristine gas (left column), and metal-rich gas (central column). Corresponding mean-squared velocity 
maps for metal-rich gas are shown on the right column. The selected region is a zoom-in on the plane z = zo ~ 370kpc//i (comoving) of the main star 
forming halo. 



3.3.1 Global features 

In Fig. [6] we show the mass profiles for metal-enriched and pris- 
tine gas of the main halo, at different redshift. The vertical lines 
refer to the virial radius at each redshift. Pollution starts at z ~ 16, 
when the early PISN spread the first metals. From this moment on, 
star formation contributes more and more to metal spreading and 
the "inside-out" mode of metal enrichment is visible in the time se- 
quence shown throughout the panels, as more and more surround- 
ing regions are affected by pollution. The structure of the pristine 
gas seems strongly influenced by star formation and feedback ef- 
fects. In fact, while star formation goes on and the metallicity of the 
environment increases, the pristine gas is always kept at ~ 0.2 kpc 
(comoving): early popIII activity halts it from getting accumulated 
toward the center. Only at redshift z < 13, when the residual popIII 
contribution to the SFR is just ~ 10 -3 , the popll regime dom- 
inates, the infall becomes more efficient and this determines the 
bump at small radii. The average amount of metal mass in the shells 
is some orders of magnitude lower than the pristine gas, but locally, 
it is possible to find highly enriched areas with metallicities up to 
Z ~ 10" 3 , i.e. - 10" 1 Z Q (e.g. see at radii of - 0.3kpc//i, in 



Fig- Winds are also able to spread around material beyond the 
virial radius, into nearby empty regions. 

This behaviour is confirmed in Fig. [7] where we show correspond- 
ing mass inflow and outflow rates for metal-rich and pristine gas, 
in the main halo. These are computed as 
dM AM 



with AM mass, Ar the radial bin, and v ra d average radial ve- 
locity for the considered shell. The pristine-gas infall rate is ~ 
10 ~ 3 — lO _1 M yr _1 . The metal outflow rates are a few orders 
of magnitude lower than the infalling rates of pristine gas, but the 
continuous ejection of metals from PISN/SN pollutes the surround- 
ings up to distances exceeding the virial radius. This is recovered 
in the inflow rates of the enriched gas. While the innermost regions 
are dominated by chaotic turbulent motions (see following discus- 
sion), the outermost enriched mass, situated at distances compa- 
rable with the virial radius, is definitely decoupled from the met- 
als in the core and is fastly moving inwards (the radial component 
of the velocity easily reaches ~ 50km/s). This incoming material 
has just been pushed away from other sources and undergone wind 
feedback. It gets into the main halo already at ~ 16, and, in the 
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following stages, polluted material systematically replenishes the 
innermost regions. Its contribution to the mass profiles becomes 
comparable or even dominant from ~ 0.5 kpc/ /i (comoving) up 
to larger radii. Since wind feedback is the main responsible for 
metal spreading up to distances larger than a few ~ kpc within 
~ 10 7 yr, the presence of metals in these dense environments af- 
ter very short timescales is a striking sign of pollution from ner- 
aby star forming haloes, whose particles are ejected and trapped in 
the neighbouring collapsed regions. The whole process, thus, leads 
to the concept of "gravitational enrichment", i.e. metal enrich- 
ment induced by different, external sources, whose polluted mate- 
rial is expelled and re-captured by the gravitational field of closeby 
haloes, loc ated at some kpc away. T his reminds of genetic enrich- 
ment (e.g. IScannapieco et al.| [2003), but, in this case, there is no 
obvious parental relationship (child/progenitor haloes) among the 
involved objects and they can have quite dinstinct evulution history. 
In contrast to self-enrichment, that is probably the main way of en- 
riching rare, isolated haloes, up to overcritical Z, the gravitational- 
enrichment mechanism appears to be the dominant one, especially 
in crowded and actively star forming environments: it can highly 
boost the metallicities and allow the gas to exceed Z cr u in only 
~ 10 7 yr (see further discussion in Sect.|4]). 
In Fig. [8j we show velocity maps for enriched and pristine gas, at 
different redshifts. The initial collapse of primordial gas appears 
along the field lines at high z (left column). This leads to star for- 
mation, followed by metal pollution, with consequent spreading of 
particles, because of wind effects. Since the typical temperatures of 
the star forming regions vary between a few hundred Kelvin (for the 
cold phase) and ~ 10 6 K (for the gas heated by SN feedback), the 
local sound speed is ~ 1 — 10 2 km/s, and, therefore, the star form- 
ing environment is characterized by supersonic flows. Moreover, 
the interplay between the hot, outflowing, enriched gas (central and 
right column), and the cold, infalling, pristine gas, leads to the de- 
velopment of a turbulent regime on kpc scales. To quantify thifl 
we give an order-of-magnitude es timate of the Reynolds number 
(e.g. see lElmegreen & Scaloll2004l . and references therein), 



R rsj ul/v rsj Miner, 



(3) 



with, u velocity, I typical scale of the flux, v kinematic viscosity of 
the gas, M Mach number, n number density, and a particle cross 
section. We have u ~ 10 2 km/s, a fiducial scale - see Figs. [6j 
and El - of I rsj kpc (in the simulations the typical dimensions of 
the star forming regions at the redshifts of interest go up to sev- 
eral comoving kpc), and v ~ c/na, where c is the sound speed. 
Typical densities of the s tar forming regions are n > 10 2 cm -3 , 
with a ~ 10 -15 cm 2 felmegreen & Scald |2004) . This implies 
v ~ 10 18 - 10 20 cm 2 /s, and R ~ 10 8 - 10 10 . Such high Reynolds 
numbers arise from the strongly turbulent state of the star forming 
regions. These are, of course, order-of-magnitude estimates, but 
they are solid enough to characterize the nature of such flows. In 
addition, the interactions between different gas phases generates 
strong hydro-instabilities: e.g., Raleigh-Taylor instabilities, origi- 
nated by the hot, enriched gas pushing against the cold, pristine 
one in its surroundings (shocks), or Kelvin-Helmholtz instabilities, 
due to the velocity differences between the outgoing, enriched gas 
and the infalling, pristine gas. The interplay among the two phases 
is very evident in the pristine-gas field lines (left column), com- 
pared with the simultaneous metal-gas ones (central column), and 

7 The effects of numerical viscosity on structure masses and velocy fields 
are wi dely discussed in li terature. See, e.g.jDolag et al.l 120050 . Vald arninil 
bOlCh . or Appendix A in lBiffi et al.1 bold) . 



with the corresponding velocity maps (right column). These con- 
siderations are expected to hold not only in this particular case, but 
also for other small primordial haloes, that are highly affected by 
feedback effects. 



3. 3. 2 PopIII and popll 

In Fig. [9j we summarize our findings about chemical feedback 
by showing in detail the radial distribution of popIII and popll 
species, in the innermost regions of the main halo, at redshift 15.00, 
13.01, and 11.00. As time passes, the most striking and evident 
feature is the gradual disappearance of the popIII generation. At 
z ~ 15, almost half of the gas is still popIII at comoving distances 
> 0.2 kpc/h. Later on, at z ~ 13, the situation is similar, but the 
residual popIII gas is found at larger and larger distances from the 
center (more than ~ 2 kpc/h comoving). Indeed, new polluted ma- 
terial has been ejected, and, as a consequence, the metallicities of 
the pre-enriched regions have increased. So, only the pristine gas 
far away from star formation sites can maintain a residual metallic- 
ity below Z cr it- By redshift ~ 11, ongoing spreading events will 
continue to confine the pristine gas to the peripheral regions, sev- 
eral kpc away from the center. Very similar trends are found when 
considering separately the infalling portion and the outflowing por- 
tion of the gas. 



4 DISCUSSION AND CONCLUSIONS 

We summarize, now, how, throughout this paper, we have addressed 
the issues presented in the introduction (Sect.Q]). We have studied 
early structure formation and metal pollution b y analyzing numer- 
ical N-body/SPH simulations, as described in llVIaio et al.l J2010h . 
With respect to o ther implementations of metal sp r eading present 
in literature (e.g . Raiteri et al. 199 61: iGnedinl 1 1998: Mosconi et al.l 
|200ll:lLia et al] 2002a b: Kawata & Gibso n 2003; Kobavashi 20041: 
Ricotti & O striker 2004; Scannapie co et al.l [2 005; Torn atore et al.l 
l200l IOpp enheimer J e7al.ll2009l : iFinlator et al.ll201ll) . our imple- 
mentation additionally follows proper molecular-chemistry treat- 
ment of e~, H, H + , H", He, He + , He ++ , H 2 , H+, D, D + , HD, 
and HeH + ; the evolution of individual metal species (C, O, Mg, S, 
Si, and Fe, since they are the most abundant element produced dur- 
ing stellar evolution); the transition from the early popIII regime to 
the subsequent popll regime; and stellar feedback for star forming 
particles. 

Metal enrichment takes place efficiently with run-away molecular 
cooling leading to the first popIII star formation events, and to the 
following, rapid transition to popll. As a consequence of the chem- 
ical feedback, higher values of Z cr u allow the popIII regime to last 
a bit longer (at most < 8 x 10 7 yr) and to pollute more with respect 
to the lower 

-Zcrit cases ( Maio et al. 2010). The first sites of metal 
pollution are clumped and dense star forming regions from which 
metals are spread away via feedback effects. 
Because of this process starting in regions with over-densities 



higher than 5 ~ 10 , the metal filling factor (Sect. 13.11 ) is ex- 
tremely small at initial times, and subsequently increases by sev- 
eral orders of magnitude, when metals get into lower-density and 
larger- volume regions (see Fig.Q] left panel). In roughly 2 x 10 8 yr, 
there is a steep increment of about 8 orders of magnitude, due to 
the explosions of the early, short-lived, massive stars and the pol- 
lution of low-density regions. Discrimination by Z cr it allows us to 
see that popIII particles dominates the filling factors only at high 
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Figure 9. Average fractions of popll (solid lines) and popIII (dashed lines) species as a function of comoving radius for the enriched material. The vertical 
dot-dashed lines indicate the virial radius. Data refer to z = 15.00 (left), z = 13.01 (center), and z = 11.00 (right). 



redshift, and are expected to be found in less dense or more iso- 
lated regions at later times. The precise value of Z cr it in the range 
[10 -6 , 10 -3 ]Z© seems to be not crucial for the properties of pri- 
mordial g as and metal en richment (Fig.Q] right panel), as also ar- 
gued in lMaio et all d201Qh . 

The phase-space evolution (Fig.[2j> and the probability distributions 
(Figs. [3] and HI) clearly indicate a dominant "inside-out" mode of 
metal pollution, according to which metals are strongly spread from 
the clumped regions to the low-density regions. Within a time of 
~ 10 8 yr, the newly enriched material accounts for ~ 90% of the 
ejected metals (see Sect. 13.21 where we compute the enrichment 
rate Ez ~ 10 _8 yr). Moreover, the fraction of enriched material 
expelled by the winds that falls in the run-away collapsing region 
and becomes available again for star formation reaches values of 
< 10% (see Sect. 13.21 ). Details about the exact value of Z cr it alter 
only slightly these findings. In order to see the effects of different 
popIII IMF's, we have also compared the top-heavy case with the 
Salpeter case. In the latter, a slightly longer delay of metal enrich- 
ment is found, due to the longer stellar lifetimes , but the subse- 
quent pollution history is quite similar (see also iMaio et alj|2010l 
for more details). 

Our dynamical analysis (e.g.: Figs. [6] |7] and [8} confirms that 
pristine-gas inflow rates of ~ 10 ~ 3 — 10 _1 M©yr _1 coexist with 
metal outflows a few orders of magnitude smaller. Metal-rich winds 
of hundreds km/s are responsible for pollution of the medium even 
beyond the virial radii, induce hydro-instabilities, and cause turbu- 
lence with Reynolds numbers of ~ 10 8 - 10 10 (Sect. [33}. 
Accumulation of metals in the high-density regions is led by "grav- 
itational enrichment" that can boost the metallicity in a few 10 7 yr 
(see e.g. Figs. [6l and [7]), without the need to wait for the growth 
of the halo to reaccrete them (self-enrichment). We note that al- 
ternative scenarios, such as self-enrichment and genetic enrich- 
ment, could not explain all the aspects shown by metal pollution 
at high redshift. For example, due to the high velocity of the ma- 
terial escaping from the star forming sites, metals can be ejected 
as far as several kpc away and the halo could effectively re- accrete 
them only if the growth process were stil l taking place. B ecause 
this requires at least some ~ 10 8 yr (e.g. iGreif et al.ll2Q10h , self- 
enrichment would fail in explaining the high m etallicities reache d 
in rsj 10 7 yr, as discussed in Sect. 13.31 (see also lMaio et alj|201oh . 
On the other hand, genetic enrichment relies on parental relation- 
ships (child/progenitor haloes), and it would not explain metal pol- 
lution in haloes which have no direct connections to each other. 
More realistically, "gravitational enrichment" can account for the 
pollution of separate haloes which evolve independently. In fact, 



when metals are spread beyond the boundaries of a halo (see e.g. 
previous sections), the gravitational force of a neighboring halo 
can attract those metals and accrete them onto it. In this way, 
metal pollution is enhanced in haloes which have gravitationally 
trapped metals, originated in completely different environments. 
The inside-out mode of metal pollution is at the origin of such dy- 
namical behaviour. 

Similar results hold for both popll and popIII components (see 
Fig. \5\ and [9]). Due to the previous considerations on metal enrich- 
ment, the popIII material is more likely to be found on the bor- 
ders of primordial haloes or in isolated, low-density environments 
(~ 50% after ~ 10 8 yr), and the popll one in more clustered re- 
gions (~ 87% after the same time). A few 10 8 yr after the onset of 
star formation, the amount of popll material becomes highly dom- 
inant over the residual popIII, reaching mass fractions higher than 
- 99%. 

One of the most striking conclusions we can draw is related to the 
rapidity and efficiency by which popll stars are formed after the 
first bursts of star formation. Because of the high metal yields of 
primordial massive PISN/SN and pollution coming from closeby 
star forming regions, the transition from the popIII to the popll 
regime is extremely rapid. Moreover, the mass trends presented in 
the previous sections (see e.g. Figs. [6l and [7} support this conclu- 
sion, as, due to star formation feedback, inflows do not efficiently 
replenish the star forming regions with newly accreted, pristine 
gas. One can estimate to first order the necessary pristine inflow 
rates, M pr ist, needed to dilute metal-rich gas, Mz, to metallicities 

Z <C Z cr it'. 

M z 



M z + M p ri. 



< Z c 
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That means (for Z -C 1): 
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M z 
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(4) 
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(6) 



i.e., for Zcrit ~ 10~ Z , one can dilute the ejected metals below 
Zcrit only with Mprist/ Mz > 5x 10 5 . These huge differences of 
several orders of magnitude between enriched-gas and pristine-gas 
rates are usually not reached, as clearly shown in Fig. [7] In general, 
this is probably what happens in the primordial haloes, which have 
similar sizes and which experience similar star formation feedback 
effects. 
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Furthermore, after the onset of star formation, the hot tempera- 
tures induced by SN thermal feedback keep dissociating molecules 
and it becomes difficult for the pristine gas to cool down: in this 
regime, the metal fine- structure transition cooling plays a domi- 
nant role, and that is the reason why one finds enriched cooling 
gas, ejected by wind feedback from nearby sources, deeper in the 
potential wells. Thus, popIII can survive mostly outside clumped 
regions, in isolated, or in low density environments, where relation 
(|6) is more easily satisfied, while popll dominates the dense, star 
forming environments. 

We note that the high resolution of our simulations and the de- 
tailed chemistry treatment allow us to make firm, consistent, and 
robust statements on the gas behaviour, from the pristine, molecular 
stages, to the following star formation phases, and the poplll/popll 
metal enrichment events. In particular, the simulations include 
a complete, self-consistent treatment of primordial gas (dealing 
with cooling and molecule evolution), which is lacking in other 
simple approximations adopted in severa l previous works (e.g . 
Kobavashi 20041: IScannapieco et al.l 120051: iTornatore et all 12007 



Cowie & McKee 1977; Brookshaw 1985; Sarazin 1988; Monaghan 


1992; Clearv & Monaghan 


1999; de Avillez & Mac Low 


2002; 


Klessen & Lin 2003; Jubelgas et al. 2004; Monaghan et al. 


2005; 



Oppenheimer et al. 2009; Finlator et al. 2011). Those implementa- 



tion can be very useful when investigating metal distributions or 
present-day star formation, but they obviously fail at high redshift, 
when molecules are the leading coolants, primordial stars form, 
and the poplll/popll transition takes place. Consistently with those 
works, metal enrichment is very patchy and inhomogeneous, with 
a characteristic "inside-out" mode. This is particularly true for the 
primordial popIII regime, whose energetics is much more powerful 
than the popll one. On the other hand, there are some caveats: given 
the small size of the boxes, we are limited in the statistical sam- 
pling of rare, large haloes. In the simulations presented here there 
is no radiative-transfer treatment from stars, which can affect the 
very close surroundings of indi vidual massive popIII stars, by heat- 
ing and ionizing th e gas (e.g. Abel et al .11 19991: lHaiman et all 2000 



Ricotti et al.1 l200ll: iKitavama & Yo shida 20 0j; llliev etaD 120051 : 
Ahn & Shapird200llwise & Abell2008l : lGreif et al.l20Q9h. but for 
low-mass SN, this effect w ould be less severe (e.g. Hasegawa et aD 
|2009|;|whalen et al.12 010). The inclusion of radiative transfer is rel- 
ativelly easy, though, only when dealing with individual or a few 
sources. It becomes very problematic and computationally expen- 
sive when studying cosmological evolution with enormous num- 
bers of stars whose radiation has to be followed simultaneously 
with gas chemistry and hydrodynamics. However, we do not expect 
substantial changes on our overall conclusions about the general 
picture of the poplll/popll transition, or the dynamics of the metals. 
Indeed, metals ejected with velocities of ~ 10 2 km/ s would prob- 
ably b e very little affected by radiative pressure (e.g. I Johnson et al.1 
120101) . The stellar radiation would probably enhance (see e.g. Fig. 
[6]) the evacuation of pristine gas from the individual min ihaloes and 
would contribute to destroy H2 molecules (e.g. lOmukai & Nishil 
1999), leaving, at the same time, the produced metals as only rele- 
vant coolants. This would essentially reproduce the same scenario 
outlined in Sect. 13.31 Star formation in purely pristine haloes de- 
layed by radiative or thermal feedback would require more mass 
to be assembled and atomic cooling to become efficient: due to the 
higher potential wells, metal spreading in such haloes would prob- 
ably be less efficient. 

Metal mixing due to diffusion processes is taken into account 
by smoothing metallicities on the kernel. This might be a c oarse 
approximation (see also discussion in iMaio et al.l 1201 (j and 
references therein), but realistic, solid and consistent imple- 
mentations of (large-scale) metal diffusion are currently lack- 
ing (approximated attempts are given by, e.g., Spitzer IS; 



Wads lev et alj bood : iGreif et al.1 120091 : Ishen et al.1 l2010h . One 
should expect that large-scale metal mixing could increase the 
metal-enriched mass which is spread around or falls back to the 
high-density environments, but no definitive conclusions can be 
drawn, at the present stage. 

The precise cosmological parameters of the ACDM model adopted 
could slightly affect the timescale of the picture, and significant 
changes in the background cosmological frame could have some 
impact, as structure formation could be signifi cantly shifted i n red- 
shift. None the less, the global trends (e.g. see Maio et al. 200^) are 
likely to be recovered. 
To conclude, we can state that: 

• metal pollution is a very efficient process and the critical 
metallicity is easily reached in less than ~ 10 8 yr; 

• the population II regime is, in general, the dominant regime of 
star formation; 

• the population III regime is expected to survive in isolated en- 
vironments or in the outskirts of star forming haloes; 

• metal pollution follows an "inside-out" mode and early metals 
are spread far away from their birth sites, even outside the virial 
radius; 

• the inside-out mode leads metal pollution in closeby star 
forming haloes, inducing "gravitational enrichment", while self- 
enrichment is a negligeble process; 

• the coexistence of cold, pristine-gas inflows and of hot, 
enriched-gas outflows determines turbulence with R ^ 10 8 — 10 10 
and hydro instabilities; 

• different Z cr it or popIII IMF alter only slightly these findings. 
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